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SIMPLIFIED  NAVIER-STOKES  EQUATIONS  AND  COMBINED  SOLUTION  OF 
NON-VISCOUS  AND  VISCOUS  BOUNDARY  LAYER  EQUATIONS 


Gao  Zhi 

(Institute  of  Mechanics,  Academia  Sinica,  Beijing) 


Abstract 


This  paper  presents  a  part  of  the  technical  report [2) 
in  which  the  author  studied  simplified  Navier-Stokes 
equations  and  the  combined  solution  of  non-viscous 
and  boundary  layer  equations.  From  the  full  Navier- 
Stokes  equations  and  an  analysis  of  the  combined  solu¬ 
tion  of  non-viscous  and  boundary  layer  equations, 
simplified  Navier-Stokes  equations  were  worked  out.  A 
perturbation  analysis  which  differs  slightly  from  the 
match-perturbation-expansions  of  inner-outer  layers 
developed  by  Van  Dyke[l]  shows  that  the  solution  of  the 
simplified  Navier-Stokes  equations  in  uniformly  valid 
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with  accuracy  of  OCRe^  )  in  the  whole  flow  field, 

where  Re^  =  is  the  density  of  free  stream, 

U^  the  x-component  of  velocity,  L  the  characteristic 
length,  the  dynamic  viscosity  of  free  steam. 


The  simplified  N-S  equations  possess  the  characteristics  of 
parabolic-hyperbolic  equations.  We  can  use  the  forward 
advancing  calculation  method  for  stationary  situations  and  it 
must  be  much  simpler  than  numerically  solving  the  elliptical 
complete  N-S  equations;  by  solving  the  simplified  N-S  equations 
we  can  simultaneously  calculate  the  non-viscous  outer  flow  and 
viscous  boundary  layer  flow.  Theoretically,  they  must  be 
superior  to  the  conventional  method  of  first  calculating  the 
non-viscous  flow  and  later  calculating  the  viscous  boundary 
layer  flow.  The  simplified  N-S  equations  can  truly  reflect  the 
mechanical  appearances  of  many  typical  flow  fields.  For  ex¬ 
ample,  they  can  accurately  calculate  the  complex  flow  chart  etc 
of  the  mutual  perturbation  between  the  viscous  boundary  layer, 
upper  entropic  layer  and  non-viscous  outer  flow  in  the  hyper¬ 
sonic  flow  field  around  a  body. 


Developments  in  the  study  of  solving  the  simplified  N-S 
equations  are  still  being  made,  in  reality,  the  problems  of 
its  mathematical  character,  steadiness,  accurate  formulation 
of  the  Cauchy  problem,  accuracy  which  can  be  reached  in  calcul¬ 
ating  the  flow  field  etc.  have  still  not  been  perfectly 
resolved  and  the  views  of  different  authors  are  not  in  agree¬ 
ment.  This  paper  presents  part  of  the  research  contents  done 
by  the  uathor  in  1967  investigating  the  simplified  N-S  equations 
brings  forth  simplified  N-S  equations  starting  from  analysis  of 
the  combined  solution  of  the  non-viscous  and  viscous  boundary 
layers  and  uses  the  perturbation  method  to  prove  that  we  can 
obtain  a  uniform  solution  with  accuracy  of  OCRe^  '  )  quantita¬ 
tive  level  from  these  equations.  It  should  be  pointed  out  in 
passing  that  analagous  analysis  for  deducing  the  simplified  N-S 
equations  were  later  seen  in  references  such  as  Reference  [3]. 

1.  Simplified  Navier-Stokes  Equations 

If  x  and  y  are  separately  the  orthogonal  coordinates  along 
the  wall  surface  and  vertical  wall  surface  (Fig.  1) ,  u  and  v 
are  the  corresponding  speed  components,  and  p  ,  p,  T ,  and  PL 
separately  indicate  the  density,  pressure,  temperature,  viscos¬ 
ity  coefficient  and  thermal  conductivity.  The  overtaking  flow 
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parameters  U*  and  fa  are  separately  the  characteristic 

quantities  of  the  speed  and  pressure  and  fixed  wall  length  L 
is  the  length's  characteristic  quantity;  we  assume  the  quantita¬ 
tive  levels  of  P  and  juu  can  use  fa  and  /fa  for  balance. 


Fig.  1 


Based  on  the  two  views  of  the  non-viscous  outer  flow  and 
viscous  boundary  layer  flow,  we  carry  out  the  following  quan¬ 
titative  estimates  of  each  term  in  the  entire  N-S  equations: 
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Key:  (1)  Non-viscous  flow;  (2)  Boundary  layer  flow; 
(3)  Non-viscous  flow;  (4)  Boundary  layer  flow. 
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Here,  n  =j  -  +  =■,  j=0  (two  dimensional  flow,  1  (axially 

r  h 

symmetric  flow) ,  H=l+ky,  k  is  the  wall  surface  curvature, 

Re^,  =  ~  '  l~  =0 ( Re-<3  1'^2)  ,  cT  is  the  boundary  layer 

thickness,  £  t  and  <£n  are  separately  the  quantitative  levels 
in  the  tangential  and  normal  equations  which  are  equal  to  or 


smaller  than  the  terms  of  the  quantita  tive 

level.  We  carry  out  the  following  discussion  when  the  large 
Reynold's  number  Re^j^l: 

1)  When  y  ,  we  omit  the  terms  with  quantitative 

levels  of  of  OiRe^  — =2 — )  in  the  entire  N-S  equations  and 
obtain  the  Euler  equations;  when  y^S*  ,  we  omit  the  terms  with 
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quantitative  levels  smaller  than  and  equal  to  OCRe^  L) 

and  obtain  the  Prandtl  boundary  layer  equations. 

2)  When  jointly  solving  the  Euler  and  boundary  layer  equa¬ 

tions  and  there  exists  mathematical  singularity  on  the  boundary 
line  of  y=<5’  ,  the  asymptotic  value  of  the  solution  of  the  Euler 
equations  when  and  the  asymptotic  value  of  the 

solution  of  the  boundary  layer  equations  when  y(  <S)  -*6* 

cannot  be  completely  the  same  and  the  singular  real  quantitative 
-l /2  n2 

1  1  •  /%  /  n  _  -i- /  \  x i.  _i  1  j  v.  _  j 3  1  »_  _  i.  _  1 .  _ 


level  is  OtRe^, 


It  should  be  noted  that  when 


y  (  <5\  the  real  quantitative  level  of  the  terms  oi 

^  —1/2 
inertia  in  the  normal  Euler  equations  is  also  0 ( Re  ' 


3)  Therefore,  in  order  to  eliminate  the  singularity  which 

appears  when  jointly  solving  the  non-viscous  outer  flow  and 

viscous  boundary  layer  flow,  we  should  select  and  omit  terms 

in  the  N-S  equations  which  are  equal  to  and  smaller  than 
-1  u^> 

°<R<io  1  )  and  afterwards  obtain  the  equations.  These  are 

then  the  simplified  N-S  equations;  by  combining  the  solution  of 
the  simplified  N-S  equations  and  Euler  equations  we  can  elim¬ 
inate  the  mathematical  singularity  and  the  obtained  accuracy  is 
the  OiRe^  )  quantitative  level  which  is  suitable  for  the 

uniform  solution  of  the  entire  flow  field. 


As  mentioned  above,  the  simplified  N-S  equations  point  to 


the  set  of  equations  obtained  after  omitting  the  terms  with 
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quantitative  levels  equal  to  and  smaller  than  0(Ret 


)  in 


the  entire  N-S  equations  estimated  according  to  the  boundary 
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layer  viewpoint.  We  concretely  write  them  (continuity  and 
energy  equations)  as: 


-Jf-  (P»r»)  +  -  (p:  Wr’>  —  0 

Ox  Oy 

/  “  0*  ,  0«  ,  $«r  ,  1  5f  .  6  /  c-.. 

'  l  w  *7  + '  aT  + -S' ;  “  -  S  67  +  s; 1 '  IT 

0y  0y  H/ 

( *  d: ■  ,  dv  ku*\  dp  ,  4  o  /  0«  \ 

'(.H67  +  ,'a7"‘F)-_aJ'  +  T6r,',ar’ 
+  J_  _9_  ('  — *  /  *2*1.  £*. 

//  0x  V  dy)  rH  dy 


0.3) 


0.0 


0.5) 
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In  the  formulas 


_  (duV  d*  n»  du  ,  .  2  0  /#*\ 

“l*)—*  Tra7  +  *“^rw 


c  is  the  constant  pressure  specific  heat.  Simplified  N-s 
P  .  _  r\  2 

s — -  term  and  when 


a.x2 


equations  (1.3)-(1.6)  do  not  include  the 

jyC^O,  it  has  four  repeated  characters  and^two  non-zero  real 

characters.  Therefore,  it  has  parabolic-hyperbolic  equation 

characteristics  [4]  and  for  constant  situations  we  can  use  the 

forward  advancing  calculation  method.  At  the  same  time,  we 

should  note  that  when  y ( <  £  )  goes  from 
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OfRe^o  )— ^O(l)  and  the  quantitative-1  levels  of  the  viscous 

terms  all  change  into  0(Re^  )  ;  accuracy  reaches  the 

0(Re  UrG>  )  quantitative  level  and  when  y^$  ,  the  simpli¬ 


fied  N-S  equations  will  smoothly  transit  to  be  Euler  equations. 
This  is  to  say  that  in  order  to  obtain  complete  flow  field 
solutions  with  the  accuracy  of  the  0(Re  J^2  )  quantitative 


•  •  • » 


level,  we  need  not  simultaneously  solve  the  simplified  N-A  and 
Euler  equations.  In  principle,  it  is  only  necessary  to  give 
the  initial  boundary  value,  which  is  the  overtaking  flow  (or 
shock  wave)  and  wall  surface  conditions,  and  independently 
solve  the  simplified  N-S  equations. 


2.  Discussion  of  Accuracy 

When  using  quasi-linear  simplified  N-S  equations  to  calcul¬ 
ate  the  entire  flow  field  including  the  non-viscous  outer  flow 
and  viscous  boundary  layer,  we  can  obtain  a  uniform  solution 
with  accuracy  of  OtRe"^^  —■ — )  .  In  order  to  prove  this  con¬ 

clusion,  it  is  necessary  to  use  slightly  different  perturbation 
expansion  from  common  inner  and  outer  lay  matched  expansion 
[1].  The  solution  of  entire  N-S  equations  (1.1)  and  (1.2), 
when  in  the  0^  y^cT  viscous  boundary  layer,  can  expand  into 
the  following  power  series  of  £ 


u/U - 

—  U„  +  £*„  +  eV,  + 

—  v,t  +  Er„  +  £JI 'a  + 

P/P- 

—  P.«  +  *P.»  +  £2p.j  + 

p/p- 

—  ft,,  +  £*,,  +  eV.a  + 

[p(X,  Y)  -  p.QOI/ep.l/1 

—  Pit  +  iPii  +  6*P.»  + 

(2.0 


In  this  series,  the  X  —  */L,  Y y/eL,  e  —  Re-*71,  —  —  0(0» 
quantities  and  their  partial  derivatives  corresponding  to  X  and 
Y  are  all  of  the  0(1)  quantitative  level.  We  substitute  expan¬ 
sion  formula  (2.1)  into  N-S  equations  (1.1)  and  (1.2),  and  by 
comparing  the  homogeneous  powers  of  £  we  obtain: 
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_  2 

In  the  formulas,  k..=kL,,  r.=r/L,  p  (X)=p  (X) /a  U  ;  we  can 
similarly  carry  out  expansion  of  the  continuity  and  energy 
equations.  It  is  necessary  to  explain  that  in  the  present 
expansion  method,  when  y  ^  (S'  ,  the  0(1)  quantitative  levels  of 
the  N-S  equations  are  close  to  being  first  order  equation  (2.2) 
and  the  0(f)  quantitative  level  equation  is  second  order 
equation  (2.3).  They  are  all  different  from  the  common  inner 
and  outer  layer  matched  expansion  method  [1] .  In  the  common 
expansion  method,  when  y  ^  (S'  ,  the  first  order  and  second  order 
normal  momentum  equations  are  separately: 

£fci.-0 

ay 

”  kjPi*1* 

Therefore,  in  the  common  expans i r.  method,  when  using  the  first 
order  equations  of  the  inner  and  outer  layers  for  simultaneous 
solution  as  well  as  when  using  second  order  equations  for 
simultaneous  solution,  there  always  exists  the  mathematical 
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singularity  of  the  O^ReS“  —  j 
boundary  line  of  y=(^. 


quantitative  level  on  the 


In  the  non-viscous  outer  flow  area,  the  solution  of 

N-S  equations  (1.1)  and  (1.2)  can  be  expanded  to  the  following 
power  series  of  £ 
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In  the  formulas,  Y  =y/L;  the  u  ,u  ,...  quantities  and  their 

e  eo  el 

partial  derivatives  corresponding  to  X  and  Y  are  all  of  the 
0(1)  quantitative  level.  We  substitute  (2.4)  into  (1.1)  and 
(1.2),  compare  the  same  order  powers  of  £  and  obtain: 
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We  can  see  that  in  the  non-viscous  flow  area  of  y  >  $  ,  the 
0(1)  quantitative  levels  of  the  N-S  equations  are  close,  that  is, 
the  first  order  equations  in  (2.5)  are  Euler  equations.  Accuracy 
reaches  the  approximation  of  the  0(£)  quantitative  level,  that 
is,  the  second  order  equations  in  (2.6)  which  are  the  modified 
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solution  after  considering  the  boundary  layer's  displacement 
thickness  effects.  From  the  combined  solution  of  equations 

(2.2)  and  (2.5),  we  can  obtain  a  uniform  solution  in  the 
entire  flow  field  with  the  accuracy  reaching  the  0(1)  quantita¬ 
tive  level;  from  the  combined  and  matched  solutions  of  equa¬ 
tions  (2.2),  (2.3)  and  (2.5),  (2.6),  we  can  obtain  the  uniform 

solution  in  the  entire  flow  field  with  the  accuracy  reaching 
1/2 

the  OfRe^  )  quantitative  level. 

By  comparing  simplified  N-S  equations  (1.4)  and  (1.5)  with 
equations  (2.2)  and  (2.3)  as  well  as  with  (2.5)  and  (2.6),  we 
can  prove  that  the  solution  of  the  simplified  N-S  equations 
satisfies  the  following  relationship  and  using  u  as  the 
example  we  have 

y  <8  I  Jj (*<•  +  *«*) j  <  0(6*  —  Re;*) 

u  (27) 

y>  8  jj C **  +  «*a)  ,  <  0(Re;‘) 

As  regards  the  other  variables,  formula  (2.7)  is  established  in 
the  same  way.  Therefore,  the  solution  of  the  simplified  N-S 
equations  is  suitable  for  the  uniform  solution  of  the  entire 
flow  field  with  accuracy  of  the  0(Re^^)  quantitative  level. 

3 .  Concluding  Remarks 

The  specific  form  of  the  simplified  N-S  equations  will  be 

slightly  different  because  of  the  differences  in  the  selection 
1/2 

of  the  OtRe^  )  quantitative  level  viscous  terms.  It  is  com¬ 
monly  considered  that  the  actual  difference  caused  by  the 

1/2 

different  selection  of  OtRe^  )  quantitative  level  viscous 
terms  is  not  large  [3-7],  However,  in  order  to  obtain  a  uni¬ 
form  solution  in  the  entire  flow  field  with  accuracy  of  the 
1/2 

0(Reoo  )  quantitative  level,  use  of  simplified  N-S  equations 

(1.3)  -(1.6)  proposed  by  the  author  [2]  is  possibly  suitable. 


The  simplified  N-S  equations  can  correctly  calculate  the 
non-viscous  flow,  high  entropic  layer  and  complex  perturbation 
between  the  boundary  layer  flows  [3-7]  in  the  shock  wave  layer. 
They  are  also  widely  used  for  the  calculation  of  other  flow 
fields:  for  example,  nonequilibrium  flows,  duct  flows,  distant 
wake  flows,  shock  wave-boundary  layer  perturbation  as  well  as 
two  dimensional  pointed  and  dull  viscous  flows  around  bodies 
etc.  The  simplified  N-S  equations  are  also  used  for  the  cal¬ 
culation  of  complex  flows  such  as  separated  flows,  hear  wake 
flows  and  those  with  compressive  corners  etc.  However,  it 
should  be  noted  that  the  simplified  N-S  equations  only  occupy 
a  dominant  position  in  viscous  effects,  that  is,  when  0, 

they  possess  parabolic-hyperbolic  equation  characteristics. 

When  in  y  viscous  effects,  they  can  be  overlooked,  that 

is,  they  transit  to  become  Euler  equations  and  the  characteristic 
root  is 


i,-i ,-H-,  a* 


>  ±  c*J  M2  —  1 
C1 


Here,  C  is  the  speed  of  sound  and  M  is  the  Mach  number.  When 
M  y  1,  .  is  the  real  character,  it  is  hyperbolic  and  the 

mathematical  type  does  not  change;  when  M<  1,  4  is  the 

repeated  character,  it  is  elliptic  and  the  type  of  the  simpli¬ 
fied  N-S  equations  changes.  The  Cauchy  problem  does  not 
pertain  to  the  simplified  N-S  equations  of  the  elliptic  type. 
It  is  essential  to  use  iterative  solutions  or  carry  out  other 
special  processing. 
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